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Abstract 

An updated analysis of all available neutrino oscillation evidence in Solar experiments 
(SK day and night spectra, global rates from Honiestake, SAGE and GALLEX) including 
the latest SNOCC and NC data is presented. Assuming that the shape of the SNO CC 
energy spectrum is undistorted and using the information provided by SNO we obtain, 
for the fraction of electron neutrinos remaining in the solar beam at energies > 5 MeV: 

^cc/^NC = 0.34t°;0t 

which is nominally 13a away from the standard value (~ 1.8 ± 0.1). The fraction of 
oscillating neutrinos which into active ones is computed to be: 

($Arc - ^cc)/i'^ssM - ^cc) = 0.92l°;39 

nearly 5a deviations from the pure sterile oscillation case. The data is still compatible 
with an important fraction of sterile component in the solar beam (up to 20% of the 
total). 

In the framework of two active neutrino oscillations we determine individual neutrino 
mixing parameters and their errors in the region of no spectrum distortion (A(Te) < 1%), 
we obtain 

Am^ = 4^.5111 X lO^^eV^ tan^ 6 = O.iOtoH. 
This is in agreement with the best solution in the LMA region. 
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1. Observation of neutral current neutrino interactions on deuterium in the SNO exper- 
iment has been recently presented Using the NC, ES and CC reactions and assuming 
the neutrino shape predicted by the SSM, the electron and active non-electron component 
of the solar flux at high energies ( > 5 MeV) are obtained. The non-electron component is 
found to be ~ 5a greater than zero, the standard prediction, thus providing the strongest 
evidence so far for flavour oscillation in the neutral lepton sector: the agreement of the total 
flux, provided by the NC measurement with the expectations implies as a by-product the 
confirmation of the validity of the SSM 

The SNO experiment measures the Solar neutrinos via the reactions |^,|^, 16]: 1) 
Charged Current (CC): + d —i- 2p + , 2) Elastic Scattering (ES): Vx + e~ ^ i^x + e~ ■ 
3) Neutral Current (NC): Ux + d p + n + Vx- The first reaction is sensitive exclusively 
to electron neutrinos. The second, the same as the one used in SuperKamiokande (SK), is 
instead sensitive, with different efficiencies, to all flavours. Finally the NC reaction is equally 
sensitive to all active neutrino species. 

The results presented recently by SNO on Solar neutrinos [jl4| confirm and are consistent 
with previous evidence from SK and other experiments g |ll|, |l^ . The CC, ES and NC 
global and day and night fluxes presented in Refs. , summarised in Table are derived 
under the assumption that the spectral shape is not distorted from the SSM prediction. 
With this assumption the SNO collaboration checks the hypothesis of non-oscillation, or zero 

(j)^+r flux. 

It would be advantageous to use this assumption of non-distortion for several reasons. 
Not only in general terms of simplicity and logical economy but also because with it a much 
higher statistical accuracy and power of prediction can be obtained (compare the total NC 
fluxes obtained with and without the distortion hypothesis obtained in |l| ) . Using these fluxes 
is appropriate for the calculation of constraints on mixing parameters in theoretical models 
where such spectrum distortion is negligible, this is true in particular for ample regions of 
oscillation space favoured by all previous data (with or without previous SNO CC data). 

Usually, the best fit to the data has been routinely obtained in the LMA region (see ||l^ 
and references therein), the qualitative explanation for that being the fact that just in that 
region the observed rather undistorted SK spectrum can be optimally adjusted. The spectrum 
distortion of the oscillating solutions for SNO in the LMA region has been explicitly found to 
be negligible [^^. Quantitavely, the main information content of the shape of the observable 
spectrum is summarised by the first moment of the distribution, the average spectrum energy. 
For the SSM case and in absence of oscillations, it is found in Ref. |37| that the average kinetic 
energy (Tg) = 7.658 it 0.006 MeV. This has to be compared with the expected value for a 
typical, non distorting, LMA oscillating solution (Tg) = 7.654 MeV and with the far values 
of the distorted SMA (7.875 MeV) and VAC (8.361 MeV) solutions. 

In this work we present an up-to-date analysis of all available Solar neutrino evidence 
including latest SNO results in the most simple framework. First we will reobtain some 
model independent results which put in a quantitative basis the extent of the deviations 
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with respect to the standard non-oscillating case and the relative importance of active/sterile 
oscillations. Second we will obtain allowed areas in parameter space in the framework of 
active two neutrino oscillations from a standard statistical analysis. Individual values for Am^ 
and tan^ with error estimation will be obtained from the analysis of marginal likelihoods. 
In this statistical analysis we include all available data from SK, Homestake and Gallium 
experiments. From SNO we include the latest results on global day and night fluxes, we 
make use in particular of: a) the total flux as measured by the NC reaction and b) the 
electron neutrino day-night global asymmetry. The main conclusion of our analysis to be 
presented below is that it is already possible to determine at present active two neutrino 
oscillation parameters with relatively good accuracy. 

2. Different quantities can be defined in order to make the evidence for disappearance 
and appearance of the neutrino flavours explicit. Letting alone the SNO data, from the three 
fluxes measured by SNO is possible to define two useful ratios, deviations of these ratios with 
respect to their standard value are powerful tests for occurrence of new physics. Here we 
compute the values for (t>cc/4>ES and (t>cc/4'NC being specially careful with the treatment 
of the correlations on the uncertitudes, the inclusion or not of these correlations can affect 
significantly their results (see table II in Ref.|Q]). For the first ratio, from the value from 
SNO rates 111 we obtain 



0CC 



0.73: 



-0.10 



a value which is ~ 2.7 a away from the no-oscillation expectation value of one. The ratio 
of CC and NC fluxes gives the fraction of electron neutrinos remaining in the solar neutrino 
beam, our value is: 

^cc Q_34+o.p5 



this value is nominally many standard deviations (~ 13a) away from the standard model 
case {(pcc/fpNC = 1.88 ± 0.08 from Ref.fS^l). 

Finally, if in addition to SNO data we consider the flux predicted by the solar stan- 
dard mode one can define, following Ref.[|l^, the quantity sin^ a, the fraction of "oscillation 
neutrinos which oscillated into active ones", again using the SNO data and fully applying 
systematic correlations (see table 2 in Ref. jl|), we find the following result: 

^SSM — S'CC 

The SSM fiux is taken as the fiux predicted in Ref.|^. Note that, although consistent 
with it, this result differs significantly from the number obtained in Ref.[18|], this is due to the 
introduction of systematic correlations in our calculation. The central value is clearly below 
one (only-active oscillations). Although electron neutrinos are still allowed to oscillate into 
sterile neutrinos the hypothesis of transitions to only sterile neutrinos is rejected at nearly 
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5(T, this significance would be reduced if we consider applying a 1-sided analysis to avoid 
non-physical values. 



3. Our determination of neutrino oscillations in Solar and Earth matter and of the ex- 



pected signal in each experiment follows the standard methods found in the literature |21| 



As is explained in detail in Ref.[12|, for this analysis we completely solve numerically the 
neutrino equations of evolution for all the oscillation parameter space. The survival proba- 
bilities for an electron neutrino, produced in the Sun, to arrive at the Earth are calculated 
in three steps. The propagation from the production point to sun's surface is computed 
numerically in all the parameter range using the electron number density rig given by the 
BPB2001 model averaging over the production point. The propagation in vacuum from 
the Sun surface to the Earth is computed analytically. The averaging over the annual varia- 
tion of the orbit is also exactly performed using simple Bessel functions. To take the Earth 
matter effects into account, we adopt a spherical model of the Earth density and chemical 
composition. In this model, the Earth is divided in eleven radial density zones in each 
of which a polynomial interpolation is used to obtain the electron density. The gluing of the 
neutrino propagation in the three different regions is performed exactly using an evolution 
operator formalism The final survival probabilities are obtained from the corresponding 
(non-pure) density matrices built from the evolution operators in each of these three regions. 
The night quantities are obtained using appropriate weights which depend on the neutrino 
impact parameter and the sagitta distance from neutrino trajectory to the Earth center, for 
each detector's geographical location. 

The expected signal in each detector is obtained by convoluting neutrino fluxes, oscillation 
probabilities, neutrino cross sections and detector energy response functions. We have used 
neutrino-electron elastic cross sections which include radiative corrections Neutrino 
cross sections on deuterium needed for the computation of the SNO measurements are taken 



from |24]. Detector effects are summarised by the respective response functions, obtained by 
taking into account both the energy resolution and the detector efficiency. The resolution 
function for SNO is that given in We obtained the energy resolution function 

for SK using the data presented in [p5|-p7|. The effective threshold efficiencies, which take 
into account the live time for each experimental period, are incorporated into our simulation 
program. They are obtained from [ |31| . 

The statistical significance of the neutrino oscillation hypothesis is tested with a standard 



X method which is explained in detail in Ref.|12]. Our present analysis is based on the 



consideration of the following quantity made of three well differentiated pieces: 

2—2,2,2 /-|\ 
X — Xglob + Xspec-sk "T XsNO- UJ 

The contribution of SNO to the is given by 
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We have introduced the flux normahsation factor a with respect to the SNO NC flux whose 
central and error values are given in Table ||. This flux normalisation will be used below as 
a scale factor for the SK spectrum. The quantity Ag is the asymmetry on the day and night 
electron neutrino rates extracted from the SNO CC, ES and NC data, imposing the condition 
Atot = ^e+/xr = 0, as predicted by active-only models 0]. The strong anticorrelation of A^ 
and Atot makes the imposition of this constraint useful, thus reducing greatly the uncertainties 
with respect to the raw Aqc- 

The definition of the Xgiob function is the following: 

x|ob = (3) 

where is the full covariance matrix made up of two terms, o"^ = cr^^^ + fcor- The diagonal 
matrix <t^^^. contains the theoretical, statistical and uncorrelated errors while a'^^j. contains 
the correlated systematic uncertainties. The _R*^'''^p are length 2 vectors containing the 
theoretical and experimental data normalised to the SSM expectations. The 2x2 correlation 



matrix have been computed using standard techniques |^, 29|. We have used data on the total 
event rates measured at Homestake experiment, at the gallium experiments SAGE fll|,|3^, 
GNO [||] and GALLEX @ (see Table (||) for an explicit list of results and references). For 
the purposes of this work it is enough to summarise all the gallium experiments in one single 
quantity by taking the weighted average of their rates. 

For the analysis of the SK energy spectrum, following closely the procedure assumed by 
the SK collaboration [^,|3^, we consider a function: 

Xs'pec = E("^*' - ^"')* (^unc + ^corCrL) ' i»R''' " R'''^ + xl • (4) 

d,n 

Where the vectors of data and expectations -Rare defined as before. We have introduced the 
SNO NC fiux normalisation factor (y, given above cuid. the correlcition pcirameter (5cor 

is assumed 

to be constrained by the last term in the sum: xj — i^cor — S^v)'^/(^s- The complete variance 
matrix is not a constant quantity. It is obtained from combining the statistical variances with 
systematic uncertainties and dependent on this correlation parameter. For each day and night 
spectrum the corresponding 19x19 block correlation matrices are conservatively constructed 
assuming full correlation among energy bins The components of the variance matrix are given 
by standard expressions [|l^] in terms of statistic errors, bin-correlated and uncorrelated 
uncertainties. The data and errors for individual energy bins for SK spectrum has been 



obtained from Ref. |31]. Other information from the SK results, such as the global day-night 
asymmetries, is to a large extent already contained in the above-mentioned quantities. It is 
therefore not included in our analysis and does not change the results presented further on. 

4. To test a particular oscillation hypothesis against the parameters of the best fit and 
obtain allowed regions in parameter space we perform a minimisation of the four dimensional 
function x^(Am^, tan^ 6*, a, (5cor)- For 5cor = S^^,0! = Omirn a given point in the oscillation 
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parameter space is allowed if the globally subtracted quantity fulfils the condition Ax^ = 
X^(Am2, 0) - x^in < xUCL). Where Xn=4(90%> 95%, ...) = 7.78, 9.4, ... are the quantiles for 
four degrees of freedom. The x^ summation now contains 41 bins in total: 3 from the global 
rates and SNO asymmetry and 2x19 bins for the SK day and night spectrums. 

The results are shown in Fig.|l](left) where we have generated acceptance contours in Am^ 
and tan^ 9. In Table (^) we present the best fit parameters or local minima obtained from 
the minimisation of the x^ function given in Eq. (^). Also shown are the values of Xmin P^^ 
degree of freedom {x^ jn) and the goodness of fit (g.o.f.) or significance level of each point 
(definition of SL as in Ref. |^^). In Fig.Q(right) we superimpose the global signal expected in 
the KamLAND experiment from reactor electron antineutrinos. We observe that at the most 
favoured regions obtained before, the KamLAND expected signal is situated at a intermediate 
region where high sensibility to both mass squared and mixing angle parameters is found. 
The experiment expects around 50% of the non-oscillation signal, in this region it would 
suffice to reach a total error of 5 — 10% (a quantity reachable after one year of data taking) 
in the measurement to be able to confirm the SNO results. 

In order to obtain concrete values for the individual oscillation parameters and estimates 
for their uncertainties, it is preferable to study the marginalized parameter constraints. It 
is justified to convert x^ likelihood using the expression C, = e~-^^/^, this normalised 
marginal likelihood is plotted in Figs. (§) for each of the oscillation parameters Am^ and 
tan^ Q. For tan^ Q we observe that the likelihood function is concentrated in a region 0.2 < 
tan^ Q < \ with a clear maximum at tan^ Q ~ 0.4. The situation for Am^ is similar. Values for 
the parameters are extracted by fitting one- or two-sided Gaussian distributions to any of the 
peaks. In the case of the angle distribution the goodness of fit of the Gaussian fit is excellent 
(g.o.f > 99.9%) even at far tail distances thus justifying the consistency of the procedure. The 
goodness of Gaussian fit to the distribution in squared mass, although somewhat smaller, is 
still good. The values for the parameters appear in Table |2[ They are fully consistent and 
very similar to the values obtained from simple x^ minimisation. 

In summary, in this work we have present an up-to-date analysis of all available Solar 
neutrino evidence including latest SNO results in the most simple framework. The direct 
measurement via NC reaction on deuterium of ^bt^ combined with the CC results has con- 
firmed the neutrino oscillation hypothesis to SOcr according to our estimate. We have obtained 
the allowed area in parameter space and individual values for Am^ and tan^ Q with error es- 
timation from the analysis of marginal likelihoods. We have shown that it is already possible 
to determine at present active two neutrino oscillation parameters with relatively good ac- 
curacy. The KamLAND and specially the near future long baseline experiments will have 
a clear chance of first, confirming present mixing parameters obtained from solar originated 
neutrinos and second, measuring first and second generation mass and mixing parameters 
under laboratory-controlled conditions. 
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Experiment [Refp | 


SssM 


Soata/SsSM (±l0') 


SK (pS8d) [30] 

n n 

SAGE [llLsa]!— 1 
GNO-GALLEX [19, 2^ 
SNO data[l,^ [| 

CC-SNO 
ES-SNO 
NC-SNO 
A,{Atot = 0) [2] 


2.32 ±0.03 ±0.08 
2.56 ± 0.16 ± 0.16 

67.2 ± 7.0 ± 3.2 
0.1 ±6.7 ±3.5 

1.76 ±0.06 ±0.09 
2.39 ±0.24 ±0.12 
5.09 ± 0.44 ± 0.45 
+0.070 ±0.052 


0.451 ±0.011 
0.332 ± 0.056 
0.521 ±0.067 
0.600 ±0.067 

0.348 ±0.020 
0.473 ± 0.053 
1.008 ±0.125 



Table 1: Summary of data used in thi^work. The observed signal {Sssm) and ratios Soata/SssM 
with respect to the BPB2001 model are reported. The SK and SNO rates are in 10^ cm~^ s~^ units. 
The CI , SAGE and GNO-GALLEX measurements are in SNU units. In this work we use the combined 
results of SAGE and GNO-GALLEX: S'ca/S'ssM (Ga = SAGE+GALLEX+GNO)=0.579± 0.050. The 
SSM «B total flux is taken from the BPB2001 model g: 0^(8B ) = 5m{lt_a^l) x 10*^ cm^^ g-i . 



Method 



A) Minimum LMA 




= 5.44 X 10-5 eV^ 


tan2 e = 0.40 xL = 30.8 


g.o.f.: 80% 


B) From Fit 


Am? 


= 4.5;2-7 X 10-5eV2 


tan2 = OAOtoll 





Table 2: Mixing parameters: A) Best fit in the LMA region and B) from fit to marginal likelihood 
distributions. 
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SNO 



CL,GA,A; 
SK Sp(d)+Sp(n) 

=7.0%+ 5.1% 
_6 <^rB)=5.09± 0.64 
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tan2 e 



KAMLAND expectations 

-2 iQ-i 1 10 

tan2 e 



Figure 1: (Left) Allowed areas in the two neutrino parameter space. The point with error bars 
corresponds to best the results from fit to marginal likelihoods. The coloured areas are the allowed 
regions at 90, 95, 99 and 99.7% CL relative to the absolute minimum. The region above the upper 
thick Hne is excluded by the reactor experiments (Right) Superimposed to the figure on the 

right, KamLAND constant signal contours normalized to the non-oscillation expectation. Contours, 
from inside to outwards, respectively at 0.4,0.5,0.6,0.7 and 0.8. 
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Figure 2: Marginalized likelihood distributions for each of the oscillation parameters Am^ (right), 
tan^ 6 (left) app The specific signature of electron antineutrinos in proton containing materials is 
the inverse beta decay process: V^+p n + e+, which produces almost isotropical monoenergetic 
positrons with a relatively high cross section. Antineutrino events would contribute to the background 
to forward-peaked neutrino solar events, earing in the fit. The curves are in arbitrary units with 
normalization to the maximum height. Dashed lines delimit ±lcr error regions around the maximum. 
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